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INTRODUCTION 

This  research  project  involves  the  analysis  of  an  existing  clinical  set  (~150  patients)  of 
frequency-domain  optical  mammograms  to  assess  the  performance  of  this  approach  to  breast 
cancer  detection.  The  analysis  of  the  breast  images  is  complemented  by  theoretical  and 
experimental  studies  to  characterize  the  proposed  algorithms  of  image  processing.  The  objective 
of  this  research  is  to  identify  the  strengths  and  the  weaknesses  of  the  current  instrument  design, 
and  to  guide  the  design  of  new  optical  instrumentation  for  breast  cancer  detection. 


BODY 

Approved  statement  of  work 

The  approved  statement  of  work  for  this  project  is  the  following: 

Task  1. 

Compute  and  analyze  the  edge-corrected  optical  mammograms  for  all  150  patients  (Months  1- 
18) 

a.  Install  the  optical  mammography  software  for  the  SUN  workstation; 

b.  Build  the  edge-effect-corrected  images  (/V-images)  from  the  amplitude  and  phase 
images; 

c.  Carry  out  a  comparison  between  the  Y-images  and  the  tumor  diagnosis  for  each 
patient; 

d.  Analyze  the  wavelength  dependence  of  the  N  parameter  and  the  ac  amplitude  for 
different  kinds  of  tumors; 

e.  Collect  a  printed  summary  of  all  the  2,400  N-images  (150  patients,  2  breasts,  2 
projections,  4  wavelengths)  and  the  spectral  features  of  all  the  tumors. 


Task  2. 

Perform  the  optical  measurements  on  breast-like  phantoms  (Months  19-24) 

a.  Prepare  the  breast-like  phantoms  (optical  inhomogeneities  +  strongly  scattering 
background); 

b.  Collect  and  analyze  the  optical  data; 

c.  Repeat  the  experiment  for  a  variety  of  sizes,  shapes,  optical  contrasts,  and 
inhomogeneous  backgrounds. 


Task  3. 

Apply  the  method  for  the  quantification  of  the  tumor  optical  properties  and  construct  the  spectral 
chart  for  benign  and  malignant  tumors  (Months  25-36) 

a.  Find  the  optical  properties,  and  the  hemoglobin-related  parameters  of  the  tumors; 

b.  Build  charts  of  benign  and  malignant  regions  in  the  ivp/  plane  at  four  wavelengths; 

c.  Determine  whether  malignant  tumors  show  a  specific  spectral  signature. 

The  first  Grant  period  (year  01)  was  devoted  to  Task  1  (which  was  planned  for  the  first  18 
months  of  the  project).  The  accomplishments  relative  to  each  one  of  the  sub-tasks  are  described 
below. 
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Task  l.a. 

The  software  provided  by  Siemens  for  the  display  of  the  optical  mammograms  has  been  installed 
on  the  SUN  workstation.  This  step  did  not  present  any  problem,  as  the  software  developed  by 
Siemens  was  compatible  with  the  new  SUN  operating  system. 

Task  l.b. 

The  software  has  been  modified  to  allow  for  the  computation  of  the  N  images  from  the 
amplitude  and  phase  images.  We  briefly  recall  that  the  so-called  N  images  are  the  result  of  an 
algorithm  of  image  processing  previously  developed  by  us  to  minimize  the  effect  of  the  breast 
thickness  variability  on  the  optical  mammograms  [Fantini  et  al.,  1996].  The  N  images  enhance 
the  tumor  detectability  with  respect  to  amplitude  or  phase  images.  The  details  of  the  image 
processing  algorithm  leading  to  the  N-images  can  be  found  in  the  appended  publication  by 
Fantini  et  al.,  1996  (Medical  Physics,  vol.23).  In  particular,  the  modified  software  implements  a 
versatile  approach  that  allows  the  user  to  select  the  values  of  the  parameters  required  to  compute 
the  N  image  (reference  pixel,  average  optical  properties  of  breast  tissue,  etc.).  A  representative 
case,  N-image  at  850  nm,  is  shown  in  Fig.  1.  The  case  in  Fig.  1  refers  to  a  71  year-old  patient 
affected  by  breast  cancer.  The  cancer,  3  cm  in  diameter,  is  clearly  detected  in  both  views  of  the 
left  breast  (LC:  left  craniocaudal;  LO:  left  oblique). 
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Fig.  1.  Representative  N-images  of  the  left  breast  of  a  71  year-old  patient  affected  by  cancer.  The  two  images  refer 
to  the  craniocaudal  (LC)  and  oblique  (LO)  projections.  The  3  cm  cancer  is  clearly  detected  in  both  views  as  a  region 
of  higher  absorbance  (dark  area). 
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The  software  displays  the  N-images  for  the  four  wavelengths  used  (690,  750,  788,  856  ran),  and 
saves  the  images  as  ASCII  files  for  further  processing.  An  example  of  the  four  N-images,  one 
per  wavelength,  of  the  same  breast  is  shown  in  Fig.  1.  This  figure  shows  that  the  cancer  (in  the 
center  of  the  image)  and  the  blood  vessels  (identified  by  their  thread-like  shape)  may  show 
different  spectral  trends.  In  fact,  in  this  case  the  cancer  contrast  decreases  with  wavelength, 
while  the  contrast  of  blood  vessels  increases  with  wavelength. 


N  left  craniocaudal 
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Fig.  2.  N  images  at  four  wavelengths  for  the  left  breast,  craniocaudal  projection,  of  an  82  year-old  patient  affected 
by  cancer.  The  cancer,  2.5  cm  in  size,  is  detected  at  the  center  of  the  image.  The  other  optical  inhomogeneities, 
darker  areas,  are  assigned  to  blood  vessels  due  to  their  morphology.  The  relative  contrast  of  the  cancer  and  the 
blood  vessels  changes  with  wavelength. 


Task  l.c. 

We  have  computed  and  printed  the  N-images  for  the  131  patients  for  which  we  had  a  complete 
pathology  report.  We  have  carefully  evaluated  these  N-images,  and  compared  the  optical 
diagnosis  with  the  actual  diagnosis.  Initially,  we  have  considered  an  optical  mammogram  to  be 
positive  only  if  an  abnormal  region  (dark  area)  appears  in  both  projections  of  the  same  breast. 
Following  this  criterion,  we  have  found  the  results  shown  in  Fig.  3. 
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Fig.  3.  Summary  of  the  positive  and  negative  optical  mammograms  analyzed  according  to  the  criterion  that  an 
optical  mammogram  is  positive  if  a  suspicious  area  is  detected  in  both  views  of  the  same  breast.  The  upper  panel 
refers  to  cancerous  breasts,  while  the  lower  panel  refers  to  healthy  breasts  and  breasts  with  benign  lesions.  The 
resulting  sensitivity  and  specificity  are  listed. 


We  have  also  applied  different  criteria,  in  order  to  build  an  ROC  (Receiver  Operating 
Characteristic)  curve  for  the  diagnostic  approach  afforded  by  single-wavelength  N-images.  The 
various  criteria  for  positivity  were  based  on  requiring  that  a  suspicious  region  be  detected  in  both 
views,  in  at  least  one  view,  with  a  contrast  beyond  a  given  threshold,  or  such  that  the  suspicious 
area  does  not  have  a  symmetrical  counterpart  in  the  contralateral  breast.  As  a  result,  we  have 
built  the  ROC  curve  shown  in  Fig.  4.  This  ROC  curve  is  very  similar  to  the  one  originally 
reported  by  the  Siemens  group  [Gotz  et  al.,  1998],  but  there  is  one  remarkable  difference.  The 
ROC  curve  of  Fig.  4  has  been  obtained  using  optical  mammograms  at  one  wavelength,  while  the 
Siemens  group  used  data  at  two  wavelengths  to  build  their  ROC  curve.  Therefore,  we  hope  to 
further  improve  upon  the  ROC  curve  of  Fig.  4,  by  exploiting  the  multi-wavelength  data  available 
to  us. 
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Fig.  4.  ROC  curve  obtained  from  the  N-images  at  one  wavelength  (690  nm)  on  a  clinical  population  of  131 
patients. 


Task  l.d. 

For  the  study  of  the  wavelength  dependence  of  suspicious  areas  detected  in  the  optical 
mammograms,  we  have  initially  implemented  a  feature  in  the  software  that  imposes  the  same 
reference  pixel  for  the  computation  of  the  N  image  at  all  wavelengths.  We  felt  that  this  is 
important  to  guarantee  that  the  data  at  each  wavelength  has  a  common  reference  point.  Then,  we 
have  developed  a  perturbation  approach,  based  on  diffusion  theory,  to  assign  an  absorption 
perturbation  to  each  optical  inhomogeneity  detected  in  the  breast.  This  perturbation  approach  is 
based  on  the  fact  that  the  absorption  perturbation  Apa  (to  within  an  unknown  wavelength 
independent  factor  C)  can  be  written  as: 


(i) 


where  ps0'  is  the  average  reduced  scattering  coefficient  of  healthy  breast  tissue,  70  is  the 
unperturbed  optical  intensity,  and  A I  is  the  measured  intensity  perturbation.  Finally,  we  have 
fitted  the  spectrum  of  the  absorption  perturbation  with  a  linear  combination  of  the  extinction 
spectra  of  oxy-  and  deoxy-hemoglobin.  By  considering  the  oxy-  and  deoxy-hemoglobin 
concentrations  ([Hb02]  and  [Hb],  respectively)  as  fitting  parameters,  we  can  translate  the 
spectral  information  into  the  tumor  oxygenation  (given  by  [Hb02]/([Hb]+[Hb02])).  This 
approach  is  depicted  in  Fig.  5  for  the  same  case  shown  in  Fig.  1. 
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Fig.  5.  The  application  of  a  perturbation  approach  assigns  the  absorption  perturbation  Apa  (to  within  an  unknown 
wavelength  independent  factor,  Q  to  a  detected  optical  inhomogeneity  (indicated  by  the  arrow  in  the  image).  By 
fitting  the  absorption  spectrum  with  the  hemoglobin  spectrum,  we  can  assign  a  tissue  saturation  value  (St02)  to  the 
detected  inhomogeneity.  The  unknown  factor  C  cancels  out  in  this  process  (because  a  multiplicative  factor  in  the 
absorption  spectrum  does  not  affect  the  tissue  saturation  value). 

We  have  applied  this  perturbation  approach  to  a  subset  of  19  patients,  obtaining  promising 
results  in  terms  of  spectral  discrimination  of  benign  and  malignant  breast  lesions.  These  results 
are  illustrated  in  Fig.  6. 
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Fig.  6.  Oxygen  saturation  of  the  hemoglobin  in  various  detected  breast  lesions  (cancers,  benign  tumors,  blood 
vessels).  The  lesion  saturation,  which  reflects  the  spectral  behavior  of  the  corresponding  absorption,  may  help 
discriminating  cancer  from  benign  lesions. 

The  result  that  malignant  lesions  have  a  higher  oxygenation  than  benign  tumors  (see  Fig.  6)  is 
surprising.  Before  trying  to  interpret  the  physiological  meaning  of  this  result,  however,  we  need 
to  better  characterize  the  perturbation  approach  used  to  obtain  Fig.  6.  To  this  aim,  we  are 
currently  performing  a  theoretical  study  and  we  are  starting  to  design  the  experimental  study 
planned  in  Task  2. 

A  difficulty  we  are  facing  is  related  to  the  fact  that  we  reference  each  N-image  to  a 
specific  pixel  in  the  mammogram  (in  fact,  a  reference  pixel  is  required  to  compute  the  N-image). 
This  means  that  the  spectral  dependence  of  N  is  not  just  affected  by  the  spectral  dependence  of 
the  tumor,  but  also  by  the  spectral  dependence  of  the  reference  pixel.  Our  first  approach  to 
tackle  this  problem  has  been  to  take  an  average  measurement  of  the  intensity  /  over  the  whole 
breast  as  a  reference  reading  (i.e.  I0  in  Eq.  (1))  for  the  perturbation  approach  described  above. 
This  may  reduce  the  variability  of  the  spectral  dependence  of  the  reference  reading.  We  are 
currently  in  the  process  of  quantifying  the  effectiveness  of  this  averaging  approach. 

Task  l.e. 

We  have  already  collected  a  printed  summary  of  the  N-images  computed  for  the  whole  clinical 
data  set  provided  by  Siemens.  Currently,  we  have  analyzed  the  spectral  information  of  a  subset 
of  19  patients,  and  we  are  currently  continuing  the  spectral  analysis,  based  on  the  perturbation 
approach,  on  a  larger  data  set.  We  are  also  considering  a  modified  perturbation  approach,  as 
well  as  other  image  processing  approaches. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Implemented  modifications  to  the  software  for  the  computation  of  N-images  and  for  the 
computation  of  the  tumor  contrast  for  each  optical  mammogram. 

•  Computation  of  the  N-images  for  the  whole  data  set  of  ~1 50  patients. 

•  Determination  of  the  ROC  diagram  based  on  one-wavelength  N  images  for  the  131 
patients  with  complete  pathology  report. 
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•  Development  of  a  perturbation  approach  for  the  spectral  analysis  of  the  optically  detected 

breast  lesions.  , 

•  Determination  of  the  effectiveness  of  the  spectral  perturbation  approach  to  discriminating 
benign  and  malignant  breast  lesions  on  a  subset  of  19  patients. 


REPORTABLE  OUTCOMES 

Manuscripts  (three  copies  of  these  manuscripts  are  enclosed  with  this  report) 

S.  Fantini,  E.  L.  Heffer,  M.  A.  Franceschini,  L.  Gotz,  A.  Heinig,  S.  Heywang-Kobrunner,  Oliver 
Schiitz,  and  Horst  Siebold,  “Optical  Mammography  with  Intensity-Modulated  Light,” 
Proceedings  Volume  from  In  Vivo  Optical  Imaging  Workshop  (September  16-17,  1999, 
National  Institutes  of  Health,  Bethesda,  MD),  A.  Gandjbakhche,  Ed.  (Optical  Society  of 
America,  Washington,  DC,  2000),  in  press. 

E.  L.  Heffer,  M.  A.  Franceschini,  S.  Heywang-Kobrunner,  L.  Gotz,  A.  Heinig,  O.  Schiitz, 
H.  Siebold,  and  S.  Fantini  “Analysis  of  Frequency-Domain  Optical  Mammograms  Using 
Spectral  Information,”  OSA  Trends  in  Optics  and  Photonics  (TOPS)  Vol.  38,  Biomedical 
Topical  Meetings,  Technical  Digest,  Postconference  Edition  (Optical  Society  of  America, 
Washington,  DC,  2000),  pp.  401-402. 


Presentations 

S.  Fantini,  “Optical  Mammography,”  Department  of  Electrical  Engineering  and  Computer 
Science,  Tufts  University,  Medford,  MA,  November  10,  1999. 

E.  L.  Heffer,  M.  A.  Franceschini,  S.  Fantini,  S.  Heywang-Kobrunner,  L.  Gotz,  A.  Heinig,  Oliver 
Schiitz,  and  Horst  Siebold,  “Analysis  of  Frequency-Domain  Optical  Mammograms  Using 
Spectral  Information,”  OSA  Biomedical  Topical  Meeting,  Miami,  FL,  April  2-5,  2000. 

S.  Fantini,  “Using  intensity-modulated  light  top  detect  breast  cancer,”  Center  for  Innovative 
Minimally  Invasive  Therapy  (OMIT),  Massachusetts  General  Hospital,  Boston,  MA,  May  2, 
2000. 

E.  L.  Heffer,  “Spectral  Analysis  of  Frequency-Domain  Optical  Mammograms,”  Beckman  Laser 
Institute,  University  of  California  at  Irvine,  May  12, 2000. 

E.  L.  Heffer,  M.  A.  Franceschini,  O.  Schiitz,  H.  Siebold,  L.  Gotz,  A.  Heinig,  S.  Heywang- 
Kobrunner,  and  S.  Fantini  “Analysis  of  Frequency-Domain  Optical  Mammograms  Using 
Spectral  Information,”  Gordon  Conference  on  Lasers  in  Medicine  and  Biology,  New  London, 
CT,  June  10-15,2000. 
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CONCLUSIONS 

We  have  completed  the  initial  analysis  of  the  optical  mammograms  based  on  the  N-images  at  the 
four  wavelengths  used.  This  analysis  has  led  to  the  ROC  curve  reported  in  Fig.  4.  We  have 
started  the  proposed  study  of  the  spectral  dependence  of  the  optically  detected  tumors.  To  this 
aim,  we  have  developed  a  perturbation  approach,  whose  application  to  a  subset  of  19  patients  has 
led  to  the  results  reported  in  Fig.  6.  We  plan  to  continue  this  analysis,  and  to  complement  it  with 
alternative  approaches. 

The  relevance  of  our  results  is  that  single  wavelength  N-images  (Fig.  4)  are  as  effective 
as  a  previously  reported  two-wavelength  analysis.  We  hypothesize  that  by  exploiting  the 
additional  information  provided  by  four-wavelength  data,  we  can  further  improve  the 
performance  of  frequency-domain  optical  mammography.  This  is  the  focus  of  our  current 
research  activity  related  to  this  project. 
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We  have  investigated  the  problem  of  edge  effects  in  laser-beam  transillumination  scanning  of  the 
human  breast.  Edge  effects  arise  from  tissue  thickness  variability  along  the  scanned  area,  and  from 
lateral  photon  losses  through  the  sides  of  the  breast.  Edge  effects  can  be  effectively  corrected  in 
frequency-domain  measurements  by  employing  a  two-step  procedure:  (1)  use  of  the  phase  infor¬ 
mation  to  calculate  an  effective  tissue  thickness  for  each  pixel  location;  (2)  application  of  the 
knowledge  of  tissue  thickness  to  calculate  an  edge-corrected  optical  image  from  the  ac  signal 
image.  The  measurements  were  conducted  with  a  light  mammography  apparatus  (LIMA)  designed 
for  feasibility  tests  in  the  clinical  environment.  Operating  in  the  frequency-domain  (110  MHz),  this 
instrument  performs  a  transillumination  optical  scan  at  two  wavelengths  (685  and  825  nm).  We 
applied  the  proposed  two-step  procedure  to  data  from  breast  phantoms  and  from  human  breasts. 
The  processed  images  provide  higher  contrast  and  detectability  in  optical  mammography  with 
respect  to  raw  data  breast  images.  ©  1996  American  Association  of  Physicists  in  Medicine . 

Key  words:  optical  mammography,  transillumination,  laser  beam  scanning,  breast  tumor, 
frequency-domain 


I.  INTRODUCTION 

Among  the  various  potential  applications  of  optical  imaging 
of  biological  tissue,  the  investigation  of  breast  cancer  is  one 
of  the  most  promising.  Breast  transillumination  with  continu¬ 
ous  wave  (CW)  light  was  first  proposed  in  1929  to  detect 
breast  lesions.1  This  CW,  broad-beam  transillumination 
methodology,  also  called  diaphanography  or  lightscanning, 
has  been  thoroughly  investigated  in  the  past  two  decades  to 
address  its  feasibility  and  effectiveness.2"8  The  main  diffi¬ 
culty  with  optical  mammography  arises  from  the  strongly 
diffusive  nature  of  light  propagation  in  breast  tissue.  This 
feature  reduces  contrast8  and  resolution9  of  optical  images, 
thus  limiting  the  sensitivity  of  the  method  for  small  and  deep 
tumors.6,8  Higher  sensitivity  can  be  achieved  with  colli¬ 
mated,  narrow-beam  transillumination,10  which  delivers  a 
narrow  interrogating  laser  beam  on  one  side  of  the  breast  and 
employs  a  localized  detector  on  the  other.11,12  Recently, 
time-resolved  methods  have  been  developed  to  provide  data 
with  higher  information  content  than  the  simple  intensity  at¬ 
tenuation  provided  by  CW  techniques.  These  time-resolved 


methods  work  in  either  the  time-domain  or  the  frequency- 
domain.  In  the  time-domain,  one  measures  the  transforma¬ 
tion  of  a  short  light  pulse  (typically  a  few  picoseconds 
FWHM)  traversing  the  tissue.13,14  In  the  frequency  domain, 
the  intensity  of  the  light  source  is  modulated  at  a  frequency 
on  the  order  of  108  Hz,  and  the  measured  quantities  are  the 
phase  shift  and  amplitude  attenuation  of  the  tissue  transmit¬ 
ted  intensity  wave.15"20  In  the  time-domain,  it  is  possible  to 
apply  a  time-gate  to  the  detector  to  select  photons  which 
have  traveled  shorter  path  lengths.21"23  Time  gating  im¬ 
proves  the  image  contrast  and  resolution.24,25  In  the 
frequency-domain,  a  similar  result  can  be  achieved  by  em¬ 
ploying  high  modulation  frequencies.26,27  However,  im¬ 
provement  in  contrast  and  resolution  is  obtained  at  the  ex¬ 
pense  of  a  decreased  signal-to-noise  ratio  due  to  the  low 
signal  from  weakly  scattered  photons.  Consequently,  to 
maximize  image  quality,  one  seeks  the  best  compromise  be¬ 
tween  the  requirements  of  improved  contrast/resolution,  and 
a  high  signal-to-noise  ratio.  In  general,  the  higher  informa¬ 
tion  content  of  time  resolved  measurements  can  be  used  to 
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Fig.  1.  Schematic  diagram  of  the  frequency-domain  light  mammography 
apparatus  (LIMA).  Al,  A2,  A3  are  radio  frequency  amplifiers;  LD1  and  LD2 
are  laser  diodes  emitting  at  685  and  825  nm,  respectively;  APD  is  an  ava¬ 
lanche  photodiode; /=  110  MHz;  A/,  =  l  kHz;  A/2=0.8  kHz. 

separate  and  quantify  absorption  and  scattering  contributions 
to  intensity  attenuation.14,28  In  the  future,  this  capability  may 
lead  to  functional  imaging,  or  to  highly  specific  optical  mam¬ 
mography.  As  an  alternative,  the  additional  information  can 
be  used  to  treat  geometrical  effects  (such  as  breast  thickness 
variation  within  the  scanned  region)  and/or  boundary  condi¬ 
tions.  In  this  paper  we  have  followed  the  latter  approach. 

This  article  presents  a  frequency-domain  light  mammog¬ 
raphy  apparatus  (LIMA)  that  performs  laser-beam  transillu¬ 
mination  of  the  breast.  Our  objective  is  to  produce  the  high¬ 
est  possible  image  contrast.  To  this  end,  we  investigate  the 
problem  of  determining  a  parameter,  derived  from  the 
frequency-domain  raw  data,  that  enhances  the  detectability 
of  optical  inhomogeneities  in  the  breast.  In  particular,  this 
parameter  should  correct  for  effects  related  to  tissue  thick¬ 
ness  variability  within  the  scanned  region,  and  to  photon 
losses  through  the  sides  of  the  breast.  We  have  performed 
laboratory  measurements  on  breast  simulating  phantoms  to 
test  the  effectiveness  of  our  derived  parameter.  Finally,  we 
have  applied  our  method  to  clinical  data  obtained  with  the 
frequency-domain  light  mammography  apparatus. 

It.  LIGHT  MAMMOGRAPHY  APPARATUS  (LIMA) 

A  schematic  diagram  of  the  frequency-domain  light  mam¬ 
mography  apparatus  (LIMA)  is  shown  in  Fig.  1.  This  instru¬ 
ment  performs  a  transillumination  optical  scan  of  the 
breast.29  The  two  laser  diodes  employed  as  light  sources  emit 
at  685  and  825  nm.  The  average  emitted  power  is  about  10 
mW.  The  laser  intensities  at  685  and  825  nm  are  sinusoidally 
modulated  at  a  frequency  of  110.0010  MHz  and  110.0008 
MHz,  respectively.  The  intensity  modulation  is  obtained  by 
direct  modulation  of  the  supplied  current  provided  by  a  DDS 
(Direct  Digital  Synthesizer)  frequency  synthesizer  and  radio 
frequency  amplifiers.  The  optical  detector  is  an  avalanche 
photodiode  (APD)  (Advanced  Photonics,  Inc.  model  SD060- 
70-62-521).  The  APD  output  current  is  sent  to  one  input  of  a 
mixer,  and  the  other  mixer  input  receives  a  signal  at  a  fre- 
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quency  of  110  MHz.  The  mixer  output  contains  low  fre¬ 
quency  components  at  1  kHz  (relative  to  the  signal  at  685 
nm)  and  at  800  Hz  (relative  to  the  signal  at  825  nm).  These 
two  low  frequency  signals  are  digitally  filtered  and  processed 
to  yield  the  frequency-domain  raw  data  (ac  amplitude  and 
phase  4>)  at  both  wavelengths.  The  separation  between  the 
signals  at  685  and  825  nm  is  obtained  electronically,  not 
optically.  The  apparatus  does  not  employ  optical  fibers.  The 
laser  beams  are  collimated  by  means  of  aspheric  and  cylin¬ 
drical  lenses,  and  dielectrically  coated  optical  prisms  make 
them  collinear.  The  APD  detector  has  a  high  numerical  ap¬ 
erture  of  approximately  0.4,  thus  providing  highly  efficient 
collection  of  transmitted  light.  In  a  typical  measurement,  the 
laser  beam  and  the  APD  detector  are  scanned  in  tandem, 
thereby  continuously  maintaining  collinearity.  The  diameter 
of  the  laser  beam  is  2  mm,  and  the  APD’s  sensitive  area  is 
1.5  mm2.  The  breast  is  compressed  between  two  transparent 
glass  plates  and  occupies  the  spatial  region  between  laser 
diodes  and  the  APD  detector  (transillumination  mode).  The 
entire  compression  assembly  with  the  two  glass  plates  can  be 
rotated  by  90  degrees  to  allow  for  data  acquisition  in  cran- 
iocaudal  and  mediolateral  projections.  The  raster  scan  is  per¬ 
formed  continuously,  line  by  line,  at  a  speed  of  about  1  cm/ 
sec.  Data  acquisition  is  unidirectional.  The  pixel  size  is 
selected  by  software.  Typically,  it  is  (l-2)mmX(l-2)mm. 
The  acquisition  time  per  data  point  is  about  100  ms,  and  the 
time  required  for  a  complete  breast  scan  is  of  the  order  of  3 
minutes. 

III.  THE  PROBLEM  OF  EDGE  EFFECTS 

A.  Statement  of  the  problem 

Transillumination  optical  scanning  produces  bidimen- 
sional  projection  images  of  the  breast.  In  pixels  close  to  the 
edge  of  this  2-D  image,  the  measurement’s  fidelity  presents 
at  least  two  problems: 

(1)  Breast  thickness  variability.  The  thickness  of  breast 
tissue,  that  coincides  with  the  separation  between  the  glass 
plates  in  the  central  part  of  the  image,  decreases  towards  the 
edge  of  the  image; 

(2)  Lateral  photon  losses.  As  the  scanner  approaches  the 
edge  of  the  breast,  additional  photons  are  lost  through  the 
sides  of  the  breast. 

These  two  problems,  that  we  collectively  designate  as 
“edge  effects,”  are  pictorially  described  in  Fig.  2.  We  found 
that  edge  effects  dominate  the  raw  data  images,  at  least  in 
our  sampling  geometry.  As  a  result,  the  presence  of  optical 
inhomogeneities,  such  as  a  tumor,  appear  in  the  optical  im¬ 
ages  as  a  slight  deformation  of  the  general  pattern  deter¬ 
mined  by  edge  effects.  Under  these  conditions,  image  con¬ 
trast  is  strongly  reduced,  and  the  sensitivity  to  tumor 
detection  is  weakened. 

B.  Approach  to  edge  correction 

Our  approach  to  the  problem  of  edge  effects  is  based  on 
two  observations.  We  dedicate  this  section  to  the  description 
and  discussion  of  these  observations. 
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Fig.  2.  Pictorial  description  of  edge  effects  in  transillumination  optical  scan¬ 
ning  of  the  female  breast.  The  geometrical  breast  thickness  along  the  source- 
detector  line  changes  during  the  scan  in  the  x-y  plane.  When  the  scanner 
approaches  the  edge  of  the  breast,  there  are  additional  photon  losses  through 
the  sides  of  the  breast.  Thickness  variability  and  lateral  photon  losses  are 
collectively  indicated  as  edge  effects. 


Observation  1.  The  measured  phase  is  linearly  dependent 
on  the  geometrical  breast  thickness  along  the  source- 
detector  line ,  and  is  not  significantly  affected  by  photon 
losses  through  the  sides  of  the  breast. 

The  implication  of  this  observation  is  that  breast  thickness 
variability  along  the  scanned  region  yields  the  dominant  con¬ 
tribution  to  the  measured  phase.  To  understand  why  this  is 
the  case,  we  should  consider  the  other  three  main  factors  that 
influence  the  phase  reading:  (a)  photon  lateral  losses;  (b) 
variability  of  the  optical  pathlength  in  air;  (c)  effects  related 
to  tissue  inhomogeneities.  A  possible  way  to  evaluate  the 
impact  of  factors  (a)  and  (b)  on  the  measured  phase  is  based 
on  the  predictions  of  diffusion  theory.  This  theoretical  model 
provides  an  adequate  description  of  light  propagation  in  bio¬ 
logical  tissue.  The  basic  diffusion  equation  is:30 

dU 

--vDV2U+vvaU=q0,  (1) 

where  U  is  the  photon  density,  v  is  the  velocity  of  light  in  the 
medium,  D  is  the  diffusion  coefficient  defined  as  1/3 p!s  (with 
fis  linear  reduced  scattering  coefficients)  (Ref.  31),  and  q0  is 
the  isotropic  source  term.  The  Green’s  function  of  Eq.(l)  for 
the  photon  flux  has  been  derived  for  uniform  media  in  a 
number  of  geometries,  in  particular  for  the  infinite  slab, 
sphere,  and  cylinder.32  Arridge  et  al  have  shown  that  the 
frequency-domain  phase  is  almost  linearly  dependent  on  the 
source-detector  separation  in  all  three  geometries.32  Further¬ 
more,  the  differences  between  the  predicted  phases  in  the 
three  geometries  are  small,  implying  that  the  additional  pho¬ 
ton  losses  in  the  spherical  and  cylindrical  geometries  do  not 
significantly  affect  the  phase.32  These  results  account  for  the 
observed  linear  dependence  of  on  tissue  thickness,  and  for 
the  insensitivity  of  the  phase  to  lateral  photon  losses  (factor 
(a)).  The  variability  of  tissue  thickness  along  the  scanned 
region  is  accompanied  by  a  variability  in  optical  pathlength 
in  air.  However,  this  latter  effect  has  little  influence  on  the 
measured  phase.  In  fact,  the  phase  shift  per  unit  distance 
traveled  in  air  (1.2°  cm-1  at  100  MHz  modulation  fre¬ 
quency)  is  typically  one  order  of  magnitude  smaller  than  the 
phase  shift  per  unit  distance  in  tissue  («*20°  cm”1  at  100 


Fig.  3.  Dependence  of  the  ac  signal  at  110  MHz  on  source-detector  separa¬ 
tion  r  for  a  medium  with  ^=0.03  cm-1  and  fi's  =  10  crrT1.  The  solid  curve 
refers  to  the  slab  geometry,  the  dashed  curve  to  the  r~a  decay  with  a- 3. 
The  two  curves  are  normalized  at  r-6  cm. 

MHz  modulation  frequency)  (factor  (b)).  More  importantly, 
the  phase  shift  in  air  is  also  linear  with  traveled  distance,  so 
that  the  overall  linear  dependence  of  phase  on  tissue  thick¬ 
ness  is  not  affected.  To  evaluate  the  importance  of  localized 
tissue  inhomogeneities,  one  should  keep  in  mind  the  order  of 
magnitude  of  a  2-3  cm  tissue  thickness  variation  effect, 
which  is  of  tens  of  degrees  in  the  phase.  The  change  in  phase 
due  to  a  strong  absorber  1  cm  in  diameter  is  of  the  order  of  1 
degree  at  modulation  frequencies  of  about  100  MHz  (Refs. 
33,34,20).  This  effect  (factor  (c))  constitutes  a  slight  pertur¬ 
bation  to  the  general  phase  pattern  determined  by  tissue 
thickness  variations. 

Observation  2.  The  edge  effects  on  the  measured  ac  am¬ 
plitude  can  be  modeled  by  a  functional  dependence  on  tissue 
thickness  (r)  as  r~a ,  with  a>  0. 

This  result  is  less  rigorous  and  general  than  that  of  obser¬ 
vation  1.  The  attenuation  of  the  ac  amplitude  with  the  thick¬ 
ness  of  an  infinite  slab  is  stronger  than  exponential.  We  have 
calculated  this  attenuation  for  a  medium  with 
fjua= 0.03  cm-1  and  /x'  =  10  cm”1  on  the  basis  of  the 
Green’s  function  reported  in  Ref.  32.  The  result  is  shown  in 
Fig.  3  for  slab  thicknesses  ranging  from  1  to  6  cm.  For  com¬ 
parison,  the  negative  power  dependence  on  thickness  (with 
a— 3)  is  also  shown  in  Fig.  3.  Lateral  photon  losses  contrib¬ 
ute  to  the  observed  reduction  in  the  increase  of  ac  signal  at 
smaller  values  of  r.  However,  there  is  also  a  contribution 
from  the  different  optical  coupling  between  tissue  and  both 
source  and  detector.  As  can  be  seen  in  Fig.  2,  the  tissue 
thickness  reduction  corresponds  to  a  change  in  the  angle  of 
light  incidence  and  to  an  increased  distance  between  tissue 
surface  and  both  source  and  detector.  This  fact,  together  with 
the  non-perfect  collimation  of  laser  beam  and  APD  detector, 
accounts  for  a  variation  in  source/tissue  and  tissue/detector 
optical  couplings.  We  have  experimentally  verified — on 
breast  simulating  phantoms — that  this  coupling  effect  is  not 
dominant  in  our  instrument.  In  our  case,  the  r~a  dependence 
of  the  ac  signal  is  thus  mainly  determined  by  lateral  photon 
losses.  Typically,  we  found  a  to  be  in  the  range  1-3.  It  is 
noteworthy  that  the  effect  of  scattering35  or  absorbing  mac- 
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Fig.  4.  Comparison  between  spatial  dependence  of  the  phase  at  110  MHz  in 
the  infinite  (dashed  line)  and  slab  (solid  line)  geometries  in  a  medium  with 
/4a=0.03  cm-1  and  /x'  =  10  cm“!.  The  spatial  coordinate  r  is  the  distance 
between  source  and  detector  in  the  infinite  geometry,  and  the  slab  thickness 
in  the  slab  geometry. 


roscopic  inhomogeneities  on  the  ac  amplitude  is  stronger 
than  that  on  the  phase.  For  a  strong  absorber  1  cm  in  diam¬ 
eter,  the  ac  reduction  is  about  20-30%  (Ref.  33,34,20). 


C.  Proposed  method  for  edge  correction 

Let  us  indicate  with  jc  and  y  the  coordinates  of  a  generic 
pixel  in  the  2-D  projection  image.  Consequently,  we  denote 
with  r(jc,y)  the  geometrical  thickness  at  pixel  (x,y).  On  the 
basis  of  observation  1  in  Section  III  B,  we  can  write: 


r(x,y)  =  r0 


&(x,y)-&s 

^max” 


(2) 


where  r0  is  the  geometrical  separation  between  the  two 
planes  of  the  compression  glass  plates  (and  hence  the  maxi¬ 
mum  thickness  of  the  breast),  $(x,y)  is  the  phase  measured 
at  pixel  (x,y),  is  the  phase  of  the  source  intensity,  and 
<£max  is  the  maximum  measured  phase  in  the  scanned  region. 
It  is  also  possible  to  recover  the  tissue  thickness  map  r(x,y) 
from  the  measured  phase  without  measuring  Qs .  This  alter¬ 
native  method  is  less  general  because  not  only  does  it  as¬ 
sume  linearity  of  d>(x,y)  on  r(jc,y),  but  it  also  assumes  a  value 
for  the  slope  d^{x,y)/dr{x,y).  The  assumed  value  for  the 
slope  is  given  by: 


d<&(.r,y) 

dr(x,y) 


(3) 


where  /  is  the  modulation  frequency  of  the  source  intensity. 
Equation  (3)  rigorously  describes  the  spatial  dependence  of 
the  phase  in  the  infinite  geometry.26,36  However,  it  is  a  good 
approximation  for  the  slab  geometry,  as  is  shown  in  Fig.  4.  In 
Fig.  4,  the  spatial  dependence  of  the  phase  in  the  slab  geom¬ 
etry  (from  Ref.  32)  essentially  coincides  with  that  in  the 
infinite  geometry  for  source-detector  distances  over  3  cm. 
The  case  reported  in  Fig.  4  assumes  typical  values  for  the 
optical  coefficients  of  the  breast  in  the  near-infrared 


Fig.  5.  Actual  (solid  line)  and  calculated  (symbols)  geometrical  thickness  of 
the  Delrin  solid  block.  The  calculation  was  carried  out  using  Eq.  (4)  at  685 
nm.  Similar  results  were  obtained  at  825  nm. 


(/za-0.03  cm-1,  /x. '  =  10  cm”1)  (Refs.  37,25),  and  a  modu¬ 
lation  frequency  of  110  MHz.  From  Eq.  (3),  the  tissue  thick¬ 
ness  in  pixel  (x,y)  can  be  written  as: 


r(x,y)  =  rQ  + 


$(x,y)-<I>I1 


I  Pa 
\2D 


1  + 


27 rf 


172 


-1 


772* 


(4) 


The  required  knowledge  of  <&s  in  Eq.  (2)  is  replaced  in  Eq. 
(4)  by  the  required  knowledge  of  the  average  optical  coeffi¬ 
cients  of  the  examined  tissue  fia  and  D. 

Once  r(x,y)  is  obtained  from  either  Eq.  (2)  or  Eq.  (4),  it  is 
possible  to  derive  a  calculated  parameter,  that  we  call  N(xyy), 
that  corrects  for  edge  effects  but  maintains  all  the  informa¬ 
tion  of  the  ac  signal,  relative  to  optical  inhomogeneities.  On 
the  basis  of  observation  2  in  Section  III  B,  we  define  N(x,y) 
as: 


N(*>y)  ra(;t,y)ac(jt,;y)  ’ 

where  ac0  is  the  ac  amplitude  measured  at  a  specific  pixel 
where  breast  thickness  is  r0,  and  ac(x,y)  is  the  amplitude 
measured  at  pixel  (x ,y).  The  value  of  the  parameter  a  de¬ 
pends  on  the  shape  of  the  examined  tissue  and  on  the  varia¬ 
tions  of  source/tissue  and  tissue/detector  optical  coupling 
within  the  scanned  region.  For  our  instrument,  phantom  data 
were  best  corrected  with  a =3  and  human  breast  data  were 
best  corrected  with  a-l.  In  this  paper,  we  present  results 
obtained  using  Eq.  (4)  to  calculate  r(x,y). 


IV.  LABORATORY  RESULTS 

As  an  initial  test,  we  verified  the  effectiveness  of  Eq.  (4) 
in  determining  the  thickness  of  a  sample  r(x,y).  We  employed 
a  solid  block  of  Delrin,  a  substance  that  shows  weak  absorp¬ 
tion  and  strong  scattering  in  the  visible/near-infrared  spectral 
region.  This  Delrin  block  has  a  central  region  6  cm  thick, 
11.7  cm  long,  and  11.4  cm  wide.  On  two  sides  of  the  block 
the  thickness  decreases,  in  a  smooth  way  on  one  side  and  by 
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Table  I.  Measured  absorption  and  reduced  scattering  coefficients  of  the 
Delrin  block. 


\  (nm) 

li a  (cm  ’) 

fi's  (cm  ') 

685 

0.019  ±  0.001 

13.2  ±  0.5 

825 

0.023  ±  0.001 

11.9  ±  0.5 

steps  on  the  other.  The  total  length  of  the  block  is  23  cm.  The 
side  view  of  the  sample  is  shown  in  Fig.  5,  where  we  also 
report  the  thickness  calculated  from  Eq.  (4).  To  apply  Eq.  (4) 
we  need  to  know  the  average  optical  coefficients  fia  and 
fi’s  of  the  sample.  We  measured  them  using  the  pre¬ 
calibrated  measurement  protocol  described  elsewhere.38  The 
results  at  the  two  wavelengths  are  reported  in  Table  I.  These 
values  of  jjia  and  fi's  are  similar  to  reported  values  for  breast 
tissue  in  vivo  in  the  same  spectral  region.37,25  The  calculated 
thickness  (from  the  phase  at  685  nm)  is  in  excellent  agree- 
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ment  with  the  actual  value  down  to  about  1  cm.  The  devia¬ 
tions  from  the  actual  yalue  are  always  less  than  3  mm.  The 
step  side  of  the  phantom  shows  a  poorer  agreement  between 
calculated  r  and  actual  r,  but  the  deviations  are  still  small.  Of 
course,  a  real  breast  does  not  present  step  variations  in  thick¬ 
ness.  A  similar  result  is  obtained  from  the  phase  at  825  nm. 
The  Delrin  phantom  also  contains  a  cylindrical  scattering 
inhomogeneity  (about  4  mm  in  diameter)  along  its  longitu¬ 
dinal  symmetry  axis.  The  edge  effects  problem  is  well  illus¬ 
trated  in  Fig.  6(a),  where  we  report  the  ac  raw  data  image  of 
the  Delrin  block  at  685  nm.  The  central  scattering  inhomo¬ 
geneity  is  invisible,  because  the  image  is  dominated  by  edge 
effects  mainly  due  to  sample  thickness  variations.  In  Fig. 
6(b)  we  show  an  optical  image  based  on  the  parameter  N 
with  a= 3.  Not  only  are  edge  effects  strongly  reduced,  but 
the  longitudinal  scattering  inhomogeneity  is  now  clearly  vis¬ 
ible.  The  pixel  size  is  2  mmX2  mm  in  both  images. 


Fig.  6.  Optical  images  of  the  Delrin  solid  block  at  685  nm.  (a)  ac  image,  (b)  N  image  with  a= 3.  The  average  optical  coefficients  of  the  sample  at  this 
wavelength  (/ua=0.019  cm-1,  /z '  =  13.2  cm-1)  match  those  of  breast  tissue.  The  sample  thickness  r  varies  along  the  v  axis.  It  increases  smoothly  from  y=0 
to  y=6  cm  (r=0  at  y=0,  r—  6  cm  at  y=6  cm),  then  it  remains  constant  (r~6  cm)  from  y— 6  cm  to  y=17.7  cm,  and  finally  decreases  by  l  cm-steps  from 
y-17.7  cm  to  y=23  cm.  Along  the  longitudinal  axis  (,r=6.4  cm)  is  located  a, cylindrical  scattering  inhomogeneity  about  4  mm  in  diameter. 
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Fig.  7.  Optical  images  at  685  nm  in  craniocaudal  projection  of  a  breast  affected  by  tumor,  (a)  ac  image,  (b)  N  image  with  a—  1. 


V.  CLINICAL  RESULTS 

To  show  the  effectiveness  of  our  approach  in  processing 
optical  data  from  human  breasts,  we  have  studied  a  number 
of  clinical  cases.  We  present  here  the  results  for  a  typical 
case,  with  the  purpose  of  pointing  out  the  kind  of  improve¬ 
ment  offered  by  N  images  with  respect  to  raw  data  ac  im¬ 
ages.  A  discussion  of  clinical  results  goes  beyond  the  pur¬ 
pose  of  this  article,  and  is  presented  in  a  separate  paper.39  To 
apply  Eq.  (4)  to  breast  data,  we  have  assumed  the  following 
average  values  of  na  and  /x '  for  breast  tissue  at  X  {  (685  nm) 
and  X2  (825  nm):  \x fl(Xj)  =0.020  cm"1,  /xa(X2)  =0.025 
cm"1,  ytt'(X1)  =  10  cm"1,  yu'(\2)=9  cm"1.  These  values 
are  consistent  with  values  reported  in  the  literature.  Fig¬ 
ures  7  and  8  report  optical  images  at  685  nm  in  craniocaudal 
and  mediolateral  projections,  respectively.  They  were  taken 


on  a  breast  affected  by  cancer  as  diagnosed  by  x-ray  mam¬ 
mography  and  histological  examination  following  needle  bi¬ 
opsy.  The  tumor  size  is  about  1.5  cm.  The  patient  is  a  33  year 
old  Caucasian  female.  Panels  (a)  of  Figs.  7  and  8  are  ac 
images,  whereas  panels  (b)  are  N  images  with  a=  1.  Pixel 
size  is  1.5  mmX  1.5  mm.  The  edge  effects  dominate  the  ac 
images,  and  the  tumor  is  not  clearly  discernible.  Image  con¬ 
trast  and  the  effectively  imaged  area  are  significantly  im¬ 
proved  in  the  N  images.  An  optical  inhomogeneity  clearly 
appears  in  the  N  images  in  the  tumor  region.  The  gray  level 
distribution  is  linear  in  both  the  ac  and  N  images,  but  only 
the  N  images  are  full  scale  (white  corresponds  to  the  smallest 
value,  black  to  the  highest).  In  the  ac  images  a  scale  reduc¬ 
tion  was  necessary  to  enhance  their  features  in  the  central 
region  of  the  breast. 
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Fig.  8.  Optical  images  at  685  nm  in  mediolateral  projection  of  the  same 
breast  shown  in  Fig.  7.  (a)  ac  image,  (b)  N  image  with  a=  1. 

VI.  DISCUSSION 

A.  Use  of  ac  signal  information 

Our  method  uses  the  phase  information  to  obtain  the  tis¬ 
sue  thickness  map  r(x,y).  The  knowledge  of  r{x,y)  is  then 
used  to  improve  contrast  and  detectability  in  the  ac  image. 
This  approach  enhances  features  present  in  the  ac  image 
without  adding  phase  information  of  optical  inhomogene¬ 
ities.  The  most  important  consequence  of  this  procedure  is 
that  the  N  image  provides  only  qualitative  information  about 
tumor  size.  The  fact  that  the  ac  image  pattern  is  insensitive  to 
the  size  of  small  optical  inhomogeneities  (below  1  cm  in 
linear  size)  was  shown  experimentally.33  For  this  reason,  our 
N  image  can  be  considered  useful  as  a  detection  tool,  but  it 
has  limited  spatial  diagnostic  significance.  As  stated  in  the 
introduction,  our  approach  is  aimed  at  improving  contrast 
and  detectability,  rather  than  resolution  and  size  evaluation. 
Another  consequence  of  using  only  ac  information  to  detect 


optical  defects  is  that  we  are  not  capable  of  recognizing  scat¬ 
tering  from  absorbing  inhomogeneites.  As  of  today,  it  is  not 
clear  whether  tumors  behave  as  stronger  absorbers  or  stron¬ 
ger  scatterers  than  healthy  tissue.  However,  a  separation  be¬ 
tween  absorption  and  scattering  effects  may  yet  result  in  an 
effective  way  of  achieving  specificity. 

B.  Use  of  single  wavelength  data 

The  parameter  N{x,y)  is  obtained  from  single  wavelength 
data.  The  detection  of  a  tumor  is  thus  based  on  the  difference 
between  its  optical  properties  and  those  of  healthy  tissue  at 
that  wavelength.  It  is  possible  to  correct  for  edge  effects  also 
by  employing  data  at  two  different  wavelengths.  In  this  case, 
the  idea  is  that  the  raw  data  (phase  and  ac)  are  influenced  by 
edge  effects  to  the  same  extent  at  the  two  wavelengths.  This 
is  a  reasonable  assumption.  By  following  this  two- 
wavelength  method,  we  were  able  to  correct  very  well  for 
edge  effects.  The  drawback  of  the  two- wavelength  approach 
is  that  the  tumor’s  detectability  is  based  on  the  different 
spectral  dependence  of  the  tumor’s  optical  properties  and 
those  of  healthy  tissue  between  those  two  wavelengths.  The 
choice  of  the  two  wavelengths  thus  becomes  critical.  It  may 
also  be  dictated  by  the  particular  kind  of  tumor  and/or  bio¬ 
logical  tissue.  The  use  of  multiple  wavelengths  could  still  be 
a  valuable  tool,  and  is  actually  being  thoroughly 
investigated.40’41  In  our  experience  to  date,  the  two- 
wavelength  method  at  685  and  825  nm  was  often  successful 
in  tumor  detection,  but  contrast  and  detectability  were  im¬ 
proved  by  the  single  wavelength  N(x,y)  map. 

C.  Significance  of  the  absolute  value  of  the  parameter 
N(x,y) 

Since  r0  and  ac0  are  constants,  contrast  and  detectability 
of  the  N  image  would  not  be  affected  by  their  absence  in  the 
definition  of  Eq.  (5).  They  simply  introduce  a  scale  factor  in 
the  N  map.  However,  we  include  them  in  the  definition  of  N 
to  study  the  information  content  of  the  absolute  value  of  N. 
In  the  absence  of  inhomogeneities,  the  N  image  ideally  as¬ 
sumes  a  constant  value  of  1.  In  reality,  the  image  is  not 
totally  flat  because  of  intrinsic  tissue  inhomogeneities  and 
imperfect  correction  of  edge  effects.  An  absorbing  or  scatter¬ 
ing  defect  causes  a  decrease  in  ac(x,y)  and  hence  a  value  of 
N(x,y)  greater  than  1.  We  found  a  good  correlation  between 
the  presence  of  tumors  and  values  of  N(x,y)  greater  than  2. 
Even  if  we  are  not  aiming  at  specificity,  it  might  be  possible 
that  a  value  of  N  higher  than  2  constitutes  a  sufficient  (but 
not  necessary)  condition  for  the  presence  of  a  tumor.  We  are 
currently  investigating  the  diagnostic  significance  of  the  ab¬ 
solute  value  of  N(x,y)  in  a  larger  number  of  clinical  cases. 

D.  Estimate  of  the  optical  coefficients  of  breast  tissue 

In  Section  III  C  we  have  presented  two  ways  to  obtain 
K*,y)  from  the  measured  phase.  The  first  one  (Eq.  (2))  re¬ 
quires  a  measurement  of  the  source  intensity  phase  <t>s .  The 
second  one  (Eq.  (4))  requires  a  measurement  or  an  estimate 
of  the  average  value  of  fia  and  /z'  of  breast  tissue.  In  this 


Medical  Physics,  Vol.  23,  No.  1,  January  1996 


156 


Fantini  et  al\  Frequency-domain  optical  mammography 


156 


article,  we  have  followed  the  second  approach  based  on  Eq. 
(4).  In  applying  Eq.  (4)  to  breast  data,  we  have  used  values 
for  fia  and  i±'s  of  breast  tissue  reported  in  the  literature.37,25 
Even  if  these  studies  report  a  small  variability  of  the  optical 
coefficients  among  individuals,  a  question  may  arise  as  to 
how  critical  is  the  choice  of  these  values?  In  other  words, 
how  much  is  the  calculated  thickness  affected  by  a  mis¬ 
judged  estimate  of  jia  and/or  /*'?  To  a  first  approximation, 
even  if  the  estimated  values  of  fia  and  /z '  are  off  by  as  much 
as  50%,  the  maximum  deviation  between  real  and  calculated 
r  is  less  than  1  cm.  However,  this  deviation  constitutes  a 
systematic  error  that  will  not  significantly  affect  the  contrast 
of  the  N(x,y)  image.  The  use  of  Eq.  (2)  will  remove  this 
source  of  error,  provided  that  the  measurement  of  can  be 
performed  with  appropriate  accuracy. 

VII.  CONCLUSION 

In  this  paper,  we  have  investigated  the  problem  of  edge 
effects  in  the  laser-beam  transillumination  scan  of  the  human 
breast.  Our  results  show  that  edge  effects  can  be  effectively 
corrected  in  the  frequency-domain,  yielding  high  contrast 
images  and  excellent  detectability.  Furthermore,  the  cor¬ 
rected  images  do  not  require  any  adjustment  of  scale  ampli¬ 
tude  or  gray  level  distribution.  The  optical  images  can  thus 
be  displayed  during  the  clinical  exam  in  real  time  with  no 
additional  manipulations.  Our  results  show  practical  feasibil¬ 
ity  of  frequency-domain  optical  mammography.  We  are  cur¬ 
rently  investigating  the  clinical  value  of  such  an  optical  tool 
by  carrying  out  a  systematic  analysis  of  optical  mammogra¬ 
phy  results  in  a  number  of  clinical  cases. 
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Abstract:  We  report  the  results  of  our  initial  analysis  of  the  frequency-domain  optical 
mammograms  collected  on  a  clinical  population  of  131  patients.  The  instrument  for 
optical  mammography  employs  four  intensity-modulated  (at  70  MHz)  laser  diodes 
emitting  at  690,  750,  788,  and  856  nm.  The  amplitude  and  phase  images  of  the  breast  are 
combined  to  produce  edge-corrected  optical  images  at  each  of  the  four  wavelengths.  The 
one-wavelength  (690  nm)  edge-corrected  images  have  been  examined  according  to  two 
different  criteria  for  the  assignment  of  positive  optical  mammograms.  Criterion  1 
requires  both  views  of  the  breast  to  be  positive,  while  criterion  2  requires  at  least  one 
view  to  be  positive.  The  resulting  (sensitivity;  specificity)  was  (72%;  52%)  (criterion  1) 
and  (88%;  30%)  (criterion  2).  We  show  how  adding  spectral  information  can  improve  on 
these  results. 

OCIS  codes:  (170.0170)  Medical  optics  and  biotechnology;  (110.0110)  Imaging  systems;  (170.4580) 
Optical  diagnostics  for  medicine;  (170.3830)  Mammography;  (1 10.7050)  Turbid  media 


1.  Introduction 

Optical  mammography  is  a  safe,  painless,  and  cost-effective  technique  that  may  become  a  useful  clinical 
tool  for  the  detection  and  diagnosis  of  breast  cancer.  It  employs  near-infrared  light  (typically  in  the 
wavelength  range  670-970  nm)  to  non-invasively  probe  the  female  breast.  Over  this  spectral  region, 
hemoglobin  is  the  dominant  intrinsic  contrast  agent.  The  hemoglobin  concentration  and  its  oxygen 
saturation  are  two  key  parameters  that  may  allow  not  only  the  detection  of  cancer,  but  also  the 
discrimination  between  malignant  and  benign  breast  lesions  using  optical  methods.  Optical 
mammography  can  be  performed  using  continuous  wave  (CW)  light  [1,2],  or  using  time-resolved 
techniques  in  the  time-domain  [3-5]  and  in  the  frequency-domain  [6-11].  Previous  implementations  of 
CW  optical  mammography  in  the  1970’s  and  1980’s  under  the  names  of  diaphanography  [12]  and 
lightscanning  [13]  failed  to  produce  a  viable  clinical  tool  because  of  their  inferior  performance 
compared  to  x-ray  mammography  [14-18],  In  the  past  decade,  the  use  of  diffusion  theory  to  model  light 
propagation  into  breast  tissue,  and  the  introduction  of  time-resolved  optical  methods  have  opened  new 
possibilities  in  the  area  of  optical  mammography.  Here,  we  present  our  approach  to  frequency-domain 
optical  mammography,  and  we  report  the  results  of  a  preliminary  analysis  of  the  optical  breast  images 
collected  on  a  clinical  population  of  13 1  patients. 


2.  Methods 


The  block  diagram  of  the  apparatus  for  frequency-domain  optical  mammography  is  shown  in  Fig.  1 . 
This  instrument  is  a  four-wavelength,  homodyne-in-quadrature  system.  The  four  wavelengths  (690, 
750,  788,  856  nm)  are  electronically  discriminated  by  modulating  the  corresponding  laser  diodes  at 
different  frequencies  ranging  from  69.50  to  70.45  MHz.  The  modulation  is  achieved  by  superimposing 
a  radio-frequency  signal  to  the  DC  current  supplies  by  means  of  bias  tees.  A  modulation  frequency  in 
the  order  of  70  MHz  is  chosen  to  achieve  a  phase  shift  per  unit  distance  of  about  207cm  in  breast  tissue. 
A  pilot  laser  at  1310  nm,  modulated  at  500  kHz,  is  used  to  probe  the  geometrical  shape  of  the  breast  and 
to  drive  protection  switches  that  prevent  the  photomultiplier  tube  detector  from  being  overexposed.  The 
emission  of  the  five  laser  diodes  is  guided  by  optical  fibers  (core  diameter:  100  pm)  to  a  fiber  coupler. 
The  output  of  the  coupler  is  an  optical  fiber  having  an  internal  diameter  of  400  pm.  The  numerical 
aperture  of  this  fiber  results  in  an  illuminated  area  larger  than  1  mm2  on  the  breast  surface.  The 
maximum  optical  power  emitted  by  the  laser  diodes  at  690,  750,  788,  856,  and  1310  nm  is  30,  10,  40, 
50,  and  5  mW,  respectively.  About  70%  of  this  optical  power  reaches  the  skin  of  the  breast  via  the  fiber 
optics.  The  maximum  irradiance  on  the  skin  is  therefore  0.7  *  135  mW/mm2.  This  irradiance  is  kept  on 
a  given  spot  by  no  longer  than  the  pixel  residence  time  (8  ms)  determined  by  the  x-y  translational  stage. 
Consequently,  the  maximum  radiant  exposure  on  the  skin  is  756  mW-ms/mm2  =  75.6  mJ/cm2,  which  is 
significantly  less  than  the  maximum  permissible  exposure  to  a  laser  beam  for  skin  (329  mJ/cm2) 
reported  by  the  American  National  Standard  for  the  safe  use  of  lasers  under  these  conditions  [19]. 
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Fig.  1.  Block  diagram  of  the  instrument  for  frequency-domain  optical  mammography.  I:  current;  f:  frequency;  A:  amplitude, 
(?:  phase;  BT;  bias  tee;  S:  shutter;  Discr.:  discriminator;  x.y  Trans.:  x.y  translational  stage:  PMT:  photomultiplier  tube;  ADC: 

analog  to  digital  converter;  MUX:  multiplexer. 


The  optical  signal  is  collected  in  a  transmission  geometry  by  an  optical  fiber  bundle  having  an 
internal  diameter  of  5  mm  (core  diameter  of  the  individual  fibers:  60  pm).  The  source  and  the  detector 
fibers  are  located  on  opposite  sides  of  the  breast.  About  90%  of  the  collected  optical  signal  is  guided  to 
a  photomultiplier  tube  detector  (for  the  690,  750,  788,  and  856  nm  signals),  and  about  10%  to  a 
photodiode  (for  the  1310  nm  pilot  signal).  The  instrumental  noise  measured  on  a  tissue-like  phantom 
without  x-y  translational  motion  and  with  an  acquisition  time  of  8  ms  is  0.15%  for  the  amplitude,  and 
0.1°  for  the  phase.  Mechanical  scanning  increases  the  instrumental  noise.  Typical  noise  levels  recorded 
in  a  human  breast  examination  (with  translational  motion,  acquisition  time  per  pixel:  8  ms,  breast 
thickness:  4-8  cm)  are  2-3%  for  the  amplitude,  and  0.5-0. 8°  for  the  phase. 

During  an  examination,  the  breast  is  slightly  compressed  between  two  glass  plates.  The  distance 
between  the  compressing  plates  is  measured  and  recorded  during  each  exam.  The  optical  image  of  the 
breast  is  acquired  by  performing  a  continuous  tandem  scan  of  the  source  and  detector  fibers  over  the  x-y 
plane  defined  by  the  compression  plates.  The  source  and  detector  fibers  remain  collinear  while  the  scan 
is  being  performed.  The  resulting  2-dimensional  projection  image  of  the  breast  is  made  of  2  x  2  mm2 
pixels.  The  breast  compression  assembly  can  be  rotated  to  allow  the  acquisition  of  breast  images  in  the 
craniocaudal  and  oblique  views,  as  routinely  done  also  in  x-ray  mammography.  The  acquisition  time 
required  to  image  the  whole  breast  is  about  2  min  per  view. 

The  amplitude  and  phase  data  are  combined  into  a  dimensionless  parameter,  that  we  called  A, 
which  is  related  to  the  optical  absorbance  of  breast  tissue  [20].  The  images  based  on  this  parameter 
(A-images)  accomplish  a  significant  correction  of  edge  effects,  which  are  determined  by  the  reduction  in 
breast  thickness  and  by  the  modified  boundary  conditions  at  pixels  close  to  the  edge  of  the  breast  [20]. 
As  a  result,  the  dynamic  range  of  the  A-images  is  significantly  less  than  that  of  the  amplitude  and  phase 
images,  and  optical  inhomogeneities  in  the  breast  appear  with  a  higher  contrast.  Furthermore,  A-images 
are  displayed  using  a  linear  gray  scale  over  the  full  range  of  data  in  the  image  (white  represents  the 
lowest  value,  black  the  highest).  This  allows  for  the  display  of  the  optical  images  of  the  breast  in  real¬ 
time,  on-line,  during  the  examination.  The  examiner  is  not  required  to  perform  any  further 
manipulations  to  the  A-images  after  acquisition. 

We  have  examined  with  frequency-domain  optical  mammography  a  clinical  population  of  131 
patients.  The  diagnosis  of  the  cases  was  determined  either  by  fine  needle  biopsy  or  by  pathology 
examination  following  surgery.  In  the  interpretation  of  the  optical  mammograms,  we  have  followed  two 
distinct  criteria.  The  first  one  (criterion  1),  previously  reported  [21],  considers  an  optical  mammogram 
to  be  positive  if  it  shows  optical  abnormalities  that  are  spatially  consistent  in  the  two  views  and  that  are 
not  thread  like.  The  second  one  (criterion  2)  requires  an  optical  mammogram  to  show  a  suspicious 
region  in  at  least  one  view  to  classify  the  case  as  a  positive  finding. 

3.  Results 

The  appearance  of  breast  cancer  in  an  optical  mammogram  is  illustrated  in  Fig.  2.  This  figure  reports 
the  craniocaudal  and  oblique  views,  at  690  nm,  of  the  left  breast  of  a  72  year  old  patient  affected  by 
invasive  ductal  carcinoma.  The  cancer  size  is  2.5  cm.  Despite  the  intrinsically  low-resolution  of  optical 
mammography,  the  optical  image  shows  this  2.5  cm-diameter  cancer  with  high  contrast.  The  optical 
contrast  is  mainly  determined  by  an  increased  hemoglobin  concentration  around  the  cancerous  tissue, 
which  significantly  increases  the  local  absorption  coefficient.  The  absorption  coefficient  of  cancerous 
breast  tissue  can  be  2-3  larger  than  the  absorption  coefficient  of  healthy  breast  tissue  [22,23].  Because 
blood  not  only  absorbs  light,  but  it  also  scatters  it,  the  smaller  scattering  contrast  observed  in  breast 
cancers  (factor  1-1.3)  [22,23]  may  also  result  from  the  increased  blood  concentration  in  the  cancerous 
area.  We  observe  that  the  A-images  are  not  able  to  discriminate  absorption  from  scattering  contrast. 
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The  dark  areas  in  the  two  images  of  Fig.  2  represent  areas  of  higher  optical  absorbance,  but  they  do  not 
provide  information  on  the  relative  absorption  and  scattering  contributions  to  the  increased  absorbance. 

We  have  examined  the  one-wavelength  (690  nm)  vV-images  from  the  clinical  population 
according  to  the  two  different  criteria  described  in  section  2.  The  results  according  to  criteria  1  and  2 
are  reported  in  Table  1 . 
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Fig.  2.  Craniocaudal  (cc)  and  oblique  (ob)  views  of  the  left  (L)  breast  of  a  72  year  old  patient  affected  by  invasive  ductal 
carcinoma.  Cancer  size  is  2.5  cm.  The  cancer  appear  with  high  contrast  in  both  views  of  the  breast. 


Table  1.  Results  of  the  analysis  of  the  690  nm  frequency-domain  optical  mammograms  from  a  clinical  population  of  13 1 

patients  according  to  criteria  1  and  2. 


Both  views 
requirement 
(Criterion  1) 

At  least  one  view 
requirement 
(Criterion  2) 

True  positive  results  (TP) 

42 

53 

True  negative  results  (TN) 

79 

52 

False  positive  results  (FP) 

73 

122 

False  negative  results  (FN) 

16 

7 

Sensitivity:  TP/(TP-HFN) 

72% 

88% 

Specificity:  TN/(TN+FP) 

52% 

30% 

Pos.  predictive  value:  TP/(TP-s-FP) 

37% 

30% 

Neg.  predictive  value:  TN/(TN+FN) 

83% 

88% 

Test  efficiency:  (TP+TN)/Total  results 

5S% 

44% 

The  results  of  Table  1  are  based  on  the  A- images  at  one  wavelength  (690  nm).  Additional  information 
can  be  obtained  by  considering  the  images  at  the  other  three  wavelengths  (750,  788,  and  856  nm).  The 
specificity  of  the  examination,  and  therefore  the  diagnosis  potential,  can  benefit  from  the  analysis  of  the 
spectral  dependence  of  the  detected  breast  lesions.  For  example,  Fig.  3  reports  the  spectral  dependences 
of  three  different  optical  inhomogeneities  observed  in  three  cases:  (1)  blood  vessel,  (2)  cancer,  and 
(3)  mastopathy.  The  corresponding  optical  mammograms  (at  690  nm,  craniocaudal  view)  are  illustrated 
on  the  right  hand  side  of  Fig.  3.  In  Fig.  3,  one  can  observe  that  spectral  features  provide  a  discriminating 
capability  that  adds  to  the  morphological  information  provided  by  single-wavelength  A-images. 
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Fig.  3.  Spectral  discrimination  of  three  structures  (a  blood  vessel,  a  malignant  tumor,  and  a  mastopathy)  observed  in  three 

optical  mammograms. 


4.  Discussion  and  Conclusion 

The  results  reported  in  Table  1  must  be  interpreted  by  considering  that  they  are  obtained  from  the 
analysis  of  one-wavelength  A-images,  which  do  not  take  full  advantage  of  the  potential  of  frequency- 
domain  optical  mammography.  The  characterization  of  the  tumor  optical  properties,  possibly  by 
reconstructing  the  absorption  and  scattering  images  of  the  breast  [9],  and  the  analysis  of  spectral 
features,  which  may  give  indications  on  the  oxygenation  of  breast  lesions  [8],  are  key  factors  in  optical 
mammography.  By  fully  exploiting  these  unique  features,  it  may  be  possible  to  improve  the 
performance  of  optical  mammography  to  a  level  that  renders  it  a  clinically  attractive  technique. 
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Abstract:  Clinical  trials  of  frequency-domain  optical  mammography  have  shown  promising  results  giving 
a  sensitivity  of  76%  and  a  specificity  of  52%.  The  sensitivity  of  this  method  is  further  investigated  by 
analyzing  all  false  negative  results.  Specificity  can  be  improved  using  spectral  information. 
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1.  Introduction 

Frequency-domain  optical  mammography  using  edge  correction  algorithms  has  produced  promising  results  for  the 
detection  of  breast  cancer  in  initial  clinical  trials.  The  reported  sensitivity  and  specificity  for  a  clinical  population  of 
130  in  these  trials  were  76%  and  52%,  respectively.  There  are  two  objectives  in  this  further  investigation  of  the 
data:  (1)  to  assess  the  sensitivity  of  the  method  more  carefully,  meaning  to  report  a  study  of  the  false  negative 
results,  and  (2)  to  improve  the  specificity  by  taking  advantage  of  spectral  measurements. 


2.  Summary 

The  prototype  for  frequency-domain  optical  mammography  is  comprised  of  four  laser  diodes  emitting  four  different 
wavelengths  (690,  750, 788,  and  856  nm),  used  as  optical  sources  that  are  modulated  at  an  RF  frequency  of  70  MHz. 
The  source  and  detector  fibers  are  located  on  opposite  sides  of  the  breast,  which  is  slightly  compressed  between  two 
parallel  glass  plates.  A  planar-tandem  scan  is  then  performed  and  the  amplitude  and  phase  information  obtained  is 
used  to  create  a  two-dimensional  projection  image  of  the  breast.  Two  views  of  the  same  breast  (craniocaudal  and 
oblique)  are  always  acquired.  Because  the  thickness  of  the  breast  is  variable,  (i.e.  it  is  thicker  in  the  center  than  at 
the  edges  between  two  parallel  plates),  an  edge  correction  algorithm,  using  the  parameter,  N,  was  used  to  improve 
image  contrast.  The  /V  parameter  uses  the  amplitude  and  phase  information  at  each  pixel,  (x,y)  as  follows: 

N(x,y)=roaco/r(x,y)ac(x,y)  (1) 

where  r0  is  the  geometrical  separation  of  the  compression  plates,  ac0  is  the  ac  amplitude  at  a  pixel  where  the  breast 
thickness  is  r0,  ac(x,y)  is  the  amplitude  measured  at  pixel  (x,y),  and  r(x,y)  is  the  breast  thickness  estimated  in  the 
selected  pixel  using  the  phase  information.  [1] 


Initial  analysis  of  these  images  for  tumor  detection  was  based  on  three  criteria  mentioned  in  [2]: 

1.  There  is  a  region  of  higher  absorption  appearing  in  the  craniocaudal  and  oblique  views  of  the  same  breast; 

2.  The  locations  of  suspicious  regions  in  craniocaudal  and  oblique  views  are  consistent  with  each  other; 

3.  The  morphology  of  the  suspicious  region  is  not  thread-like  in  at  least  one  view  of  the  breast.  The  thread-like 
appearance  is  indicative  of  blood  vessels,  which  would  absorb  light  at  the  wavelengths  being  used. 
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The  results  of  this  analysis  were  then  compared  with  x-ray  mammography  data,  as  well  as  biopsy 
information  and  pathology  reports  in  order  to  determine  which  results  were  false  positive,  false  negative,  true 
positive  or  true  negative.  In  order  to  improve  the  sensitivity  of  the  optical  method,  it  is  necessary  to  reduce  the 
number  of  false  negative  results,  and  in  order  to  improve  the  specificity  it  is  necessary  to  reduce  the  number  of  false 
positive  results. 


It  has  been  determined  that  false  results  may  have  several  different  causes.  furst,  some  of  the  images  to  be 
analyzed  showed  common  patterns  that  likely  resulted  from  geometrical  effects.  This  geometrical  pattern  made  it 
difficult  to  determine  whether  there  was  actually  an  absorbing  body  within  the  breast  or  not.  In  order  to  correct  for 
this,  we  introduce  a  correction  factor  for  the  inter-subject  variability  of  the  ac  dependence  on  thickness.  This 
modifies  the  value  for  the  N  parameter,  resulting  in  //-images  with  less  evident  geometrical  patterns. 


Second,  there  were  several  cases  in  which  a  tumor  appeared  in  one  view,  but  not  in  the  other.  The  reason 
for  this  is  not  yet  known,  but  we  believe  that  it  is  related  to  the  tumor  depth  and  to  geometric  effects.  Another  cause 
of  false  results  is  the  third  criterion  listed  above,  whose  analysis  is  subjective.  An  objective  identification  of  blood 
vessels  in  the  optical  mammograms  would  reduce  the  subjectivity. 


In  several  of  the  cases  where  a  false  negative  result  was  previously  reported,  the  tumor  was  located  very 
close  to  the  chest  wall  and  it  was  outside  of  the  field  of  view. 


One  method  of  reducing  the  number  of  false  positive  results  is  to  use  spectral  information  in  order  to 
discriminate  among  blood  vessels,  benign  lesions,  and  malignant  tumors.  It  has  been  observed  that  the  absorption 
properties  of  blood  vessels  differ  from  those  of  tumors.  While  it  is  yet  to  be  determined  whether  there  is  a  consistent 
difference  in  the  spectral  properties  of  malignancies  and  benign  tumors,  encouraging  results  in  this  direction  have 
been  reported.[3] 
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